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2-D  Dynamic  Cavity  Expansion  Model 
for  Arbitrary  Tractions 


Hyung  Je  Woo 


Abstract 

A  two  dimensional  dynamic  cavity  expansion  model  is  developed  for  the  cylindrical 
cavity  in  an  infinite  medium  subjected  to  an  asymmetric  constant  velocity  vector  as  a 
surface  traction  on  the  cavity  wall.  From  the  theory  of  elasticity,  the  equation  of  motion  for 
an  elemental  volume  of  elastic  material  in  plane  strain  and  in  polar  coordinates  is  used.  The 
analysis  dealing  with  the  exterior  of  the  cavity  consists  of  Hankel  functions  of  the  first 
kind.  Two  different  methods,  the  Least  square  procedures  and  the  Fourier  transforms,  are 
employed  for  the  elastic  problems.  As  a  first  technique,  Least  square  procedures  has  been 
used  to  solve  the  unknown  constants  of  the  elastic  solutions  with  the  initial  conditions  on 
the  cavity  surface.  The  frequency  equations  can  be  solved,  for  example,  for  the  first  ten 
roots  for  each  of  the  first  eleven  modes  of  vibration.  Secondly,  Fourier  transforms  are 
used  to  express  all  the  solutions  in  the  frequency  domain  and  then  the  complete  solutions 
can  be  obtained  by  inversing  these  Fourier  transforms.  Algorithms  of  inverse  fast  Fourier 
transforms  (FFT)  will  be  an  efficient  technique  for  this  case.  For  each  discrete  time  step, 
the  equilibrium  equations,  along  with  the  von-Mises  yield  condition,  are  solved  for  the 
solutions  in  the  plastic  region.  The  model  can  be  also  applied  when  arbitrary  normal  and 
tangential  tractions  acting  on  the  surface  of  the  circular  cavity  are  prescribed. 


1.0  Elastic  Region 

1.1  Wave  Equations 

The  equations  of  motion  of  an  isotropic  elastic  solid  in  which  body  forces  are 
absent,  are: 


d2u 

d2v 

V 


d 2 


w 


dt 2 


=  (A+/r)-^-+jUV2w 
dx 

(la) 

=  (A+/r)^+/rV2v 
dy 

(lb) 

=  (X  +  H)^  +  /N2w, 
dz 

(lc) 

d2  d 2  d 2 

where  the  operator  V2  is  written  for  (— r  +  -r— r  +  tt) 

dxz  dy  dz 


Using  the  vector  operators  with  the  displacement  u,  we  have  another  form  of  the  vector 
displacement  equation  of  motion  as  shown  below: 


1 


(A  +  2  p)V(V  •  u)  -  pV  x  V  x  u  =  pii 


(2) 


where  the  dilatation  is  A  =  V  u. 

These  equations  of  motion  may  be  shown  to  correspond  to  the  propagation  of  two  types  of 
waves  through  the  medium 


P~^T  =  (A  +  2p)V2A 

at 


(3a) 


p|?  =  pV2w,  (3b) 

where  A  is  the  dilatation  propagating  through  the  medium  with  the  velocity  [(A  +  2 p)  /  p]72 

and  w  is  the  rotation  about  the  axis  propagating  with  the  velocity  (p  /  py2 .  In  vector 
forms, 


V2A  =  A  tc2  (4a) 

V2w  =  w  /c2.  (4b) 

Hence  the  displacement  equation  of  motion  governs  waves  of  the  dilatation  and  the 
rotation. 

Poisson  gave  a  class  of  particular  solutions  to  Eqn.  (2),  based  on  the  assumption 
that  the  displacement  vector  u  is  the  gradient  of  a  potential  function,  i.e.,  u  =  V(|> .  Later, 
Lame  gave  the  general  solution  using  both  a  scalar  and  vector  potential,  such  as  every 
displacement  vector  field  of  the  form 


u  =  V0  + Vx  i/A,  (5) 

which  satisfies  Eqn.  (2),  provided  p  and  yf  are  solutions  of: 

V20  =  0  /  c2  (6a) 

V2yr=  ylc2.  (6b) 

When  solutions  are  found  to  Eqn.  (6)  for  (j)  and  yr,  displacement  vector  u  can  be  obtained 
from  Eqn.  (5),  and  this  is  a  solution  to  the  displacement  equations  of  motion. 

We  particularly  consider  the  2-dimensional  plane  strain  problem  and  use  the 
cylindrical  coordinates  (r,  9,z ).  In  the  absence  of  the  body  force,  two  scalar  wave  equations 
are: 


V20  =  <p  1  cd2 

(7a) 

V>=  ifr/c;. 

(7b) 

In  cylindrical  coordinates, 
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For  harmonic  waves  of  frequency  to,  we  write 


<j)(r,e,t)  =  0(r,d)e-ia'  (8a) 

V v{r,e,t)='¥{r,0)e-i°n.  (8b) 

Since  the  potentials  must  be  periodic  in  6  with  period  2k,  we  seek  solutions  of  the  form: 

<P(r,6)  =  R(r)e‘ne  (9a) 

x¥{r,e)  =  S{r)eine  (9b) 


where  n  =  integer. 

Substituting  Eqn.  (8)  into  Eqn.  (7),  we  obtain  two  Bessel  equations: 


R"  +  -R'  +  (a2-^T)R  =  0 
r  r 

(10a) 

S"  +  -S'  +  (P2~^-)S  =  0 
r  r 

(10b) 

where  a  =  —  and  p  =  — . 

Cd  Cs 

The  general  solutions  of  these  equations  are  the  linear  combinations  of  Bessel  functions  Jn 
and  Yn,  e.g.. 

R(r)  =  AJn  (ar)  +  BYn(ar) 

(ID 

where  the  Bessel  functions  of  the  first  and  second  kind  and  order  n  are  given  by  (when 
p  =  n  =  0  or  integer): 

,  rr,  f  (-l)m(*/2)2m+p 
h  m!(m  +  p)l 

(12a) 

^w=i<-1,IW2>7 

(12b) 

,  Jp(x)cos(pjr)-J_  (x) 

Yp(x)  =  .  ,  .  P  • 

sin  (p7t) 

(12c) 

For  brevity  we  employ  the  notation: 


(13) 


C  (z)  =  J„(z) 

C?\z)  =  Yn(z ) 

C?(z)  =  H?\z)  =  Jn(z)  +  iYnW 

C?(z)  =  H?(z)  =  Jn(z)~iYn(z) 

where  Hf(z)  and  H™(z)  are  the  Hankel  functions  of  the  first  and  second  kind, 
respectively.  The  functions  Jn  and  Yn  are  suitable  for  the  solutions  of  the  interior  boundary 
value  problems,  while  Hf  and  H'2)  are  suitable  for  the  exterior  boundary  value  problems. 

When  a  cylindrical  cavity  is  subjected  to  dynamic  surface  tractions  along  the  axis  of 
the  cylinder,  the  state  of  strain  is  planar.  The  general  solution  of  an  arbitrary  cross  sectional 
cavity  problem  in  an  infinite  medium  will  have  the  Hankel  function  Hf  for  the  cylinder 
function  that  is  to  be  retained  in  the  expressions  of  the  Lame  potentials. 


Figure  1.  Cylindrical  cavity  with  surface  tractions. 


The  use  of  Hm  ensures  that  the  Sommerfeld  radiation  condition  is  satisfied  at  infinity,  so 

n 

that  H"\r)exp(-icot)  represents  waves  advancing  toward  infinity  from  the  surface  of  the 
cavity.  Note  that  for  waves  converging  to  the  cavity,  one  must  employ  Hf\r)cxp(—icot). 
Thus  for  the  problem  of  a  cavity  subjected  to  surface  loads  we  use  Cf  =  if® . 
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(14a) 


0(r,d)  =  ^AH{:\ar)eM 

n= 0 


't'(r,e)  =  YJBM\Me>ne 

/i=0 

(14b) 

Generally,  we  have 

(15) 

Vx*'-<7  » 

(16) 

£ 

and  since  u.  =  0  and  — 
dz 

=  0, 

1  dy^ 
r  dr  r  de 

(17a) 

_  1  d</>  dy/z 

6  r  dd  dr 

(17b) 

The  displacement  field,  aside  from  a  factor  exp(-tcat),  becomes  (see  Appendix  A): 

«,=-££  lAM3)(a,lr)+iB^\e,lr)]e“  (18a) 

r  «=0  k= 1 

%  =-££[M,,V;3>(a^r)  +  B„V<2I>(fcr)]e"  (18b) 

r  n=0  k= I 

where 


Uf\ar)  -  arCH\(ar) -  nCf{ar) 

Uf{M  =  +nCf(pr) 

Vl‘\ar)  =  +nC‘0(ar) 

V®  (fir)  =  -j8rC®(j8r)  +  nC^ifir). 


(19) 


In  deriving  Eqn.  (19)  we  have  employed  the  recurrence  relations  for  the  cylinder  functions 
Cn(z),  where  z  is  complex  and  n  is  any  number  (not  necessarily  an  integer): 
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zC'(z)  =  nCn(z)  -  zCn+l(z )  =  -nC„(z)  +  zC^z) . 


(20) 


The  physical  components  of  the  stress  tensor  follow  from  Hooke’s  law,  namely, 

,d2<p 


cjr  =  AV  + 


dr2 


1  dy_  1  dy_ 

r 2  d0  r  drdO 


_  ^  ^  d(j>  l  d2(l>  1  dw__d^_ 

<T,-av^+— {-+-—+-3e  Me> 


Gz  =  AV20 

.2  <?20  2  d(j)  1  <?V  ,  1 

"-■rt7*55“7a5+7a?”p-+;ir) 


(21) 


°*=0 

E  ,  ,  v£  2  v 

where  u  = - and  A  =  — - — — — —  =  - — —  [i- 

2(1+ v)  (1  +  v)(l  — 2v)  1  — 2  v 

Thus,  aside  from  the  exponential  factor  exp(-ituf),  we  have  (Appendix  B): 

r  n=0  k= I 

a.  =  #££[4*0<V)+ 

r  n=0 *=1 

r  «=0*=I 

ff*  =  2fXEK»0<V>  +  (A.Ok" 

/t=0  fc=l 

CT^=0 

where 


T(x[\ar)  =  (n2  +  n  - V  )C%\ar)  -  arC^ ( ar ) 


(22) 
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TV(M  =  +[-(n2  +  n)C{^  ((3r )  +  nfirC„\ {fir )] 

Twicer )  =  ~{n2  +n-  —fi2r2  - a2r2)C^{ar)  +  arC^iar) 

2  (23) 

^OSr)  =  +[(n2  +  n)C:°(/3r)  -  «/?<!>  (j8r)] 

^(ar)  =  +[-(n2  +  n)C^(ar)  +  narC^(ar)] 

J^08r)  =  -(n2  +  n  -  |j32r2  )C«\fir)  +  firC <«  (0r) 

7^'1)(ar)  =  (a2r2  -^•j82r2)C^,)(ar). 


1.2  Frequency  Equations 

Boundary  conditions  for  free  vibrations  are: 

crr  =  crffl  =  0  on  r  =  a.  (24) 

We  have 

+  (25a) 

a  n=0*=l 

o„  =  r<:>(arta)  +  (25b) 

a  n=0*=l 

Setting  <7r  and  <Jre  equal  to  zero  at  r  =  a,  we  obtain  a  set  of  two  homogeneous  equations. 
The  determinant  of  this  set  must  vanish  for  the  existence  of  a  non-trivial  solution  for  Ank 
and  Bnk.  Thus  we  have  the  frequency  equation: 


A  =  T^(a’*a>  (26) 

TS'(a^l)Tg>(ff^a)  +  r*\anta)Tgtfnta)  =  0  (27) 


where 

T\\\aa)  =  {n2+n- 1 fi2a2)H^\aa )  -  OtaH"\{aa) 

Tn  (Pa)  =  ~(n2+  n)H^(fia)  +  nfiaH™  0a) 
Tf]\aa)  =  -(n2  +n)Hl'\aa)  +  naaH^l(oca) 


(28) 
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72>0 to)  =  -(n!  +  "  -  ^fa2)H'"0a)  +  paH^,(fia). 


1.2.1.  Axially  Symmetric  Case:  n  =  0 

We  have: 


Tlf\aa)  =  -■ij32a2Ho1)(«a)  +  aaH[X){aa) 


T^ipa)  =  0 

T™(aa)  =  Q> 


1 


7^  (/fa)  =  jp2a2H"(l3a)  -  j8afl?,}0fc). 


Then  the  frequency  equation  becomes 


Tt?(,a,ta)T$(J3^a)  =  0. 


(30) 


For  the  given  Poisson’s  ratio,  v,  frequency  ratio  k  is  fixed,  i.e.: 


and  so, 


a2  _  fi  _  l-2v 
/F"  ~  ~X  +  2fl  ~  2(1  -  v) 


a  = 


|  l-2v 

pa- v) 


■p. 


(31) 


(32) 


If  v  =  — ,  then  k  =  — . 

4  3 

i.e.  (5  =  a/3 a 

Then  the  frequency  equation  becomes: 

[tfpo^a)-— HPP0*a)][//pV3a0ta)  — ^ — tf,(1)(V3a0,a)]  =  0.  (33) 

3  aota 

Eqn.  (33)  can  be  solved  for  the  roots  aot,  which,  in  turn,  yield  ft)0*,  if  a,  v  and  c,  are 
known. 
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1.2.2.  Nonaxially  Symmetric  Case:_n  ^  0 


The  frequency  equation  becomes,  if  v  =  — , 

TS\anta)T$(j3aMa) + T‘?(a„,a)7«)(V3a,lia)  =  0  (34) 

where,  n  =  1,2,...,N  and  k  =  0,1,2,. ..,K. 

We  have  K  number  of  frequency,  a„k,  for  each  n. 


1.3  Least  Square  Method 

Unknown  constants  Ank  and  Bnk  can  be  determined  from  the  initial  conditions  at 
t  =  0.  For  the  constant  velocity  initial  conditions. 


«r(c,0,0)  =  Vr  =  F0 


Me(fl,0,0)  =  Ve  =  0. 


(35) 


y 


Figure  2.  Velocity  initial  condition. 


From  Eqns.  (18),  radial  and  circumferential  velocities  are  the  real  part  of: 

4",  =  “(«.»'•)-  (36a) 

r  n= 0 k= 1 

4“.  =1XX[Aa®„Vl<J,(a„r)-«rt<BJlV<3»(/3 ,tr)kMe-'“  (36b) 

01  r  n= 0  *=1 
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Then 


Re  V,  =  -  X  Z  [A,*  (a nkr)  sin  n6  -  BllkconkU^(Jinkr)cosne]  (37a) 

r  n-0  k= 1 

ReVe  =  ^SS[4,^^(n^)cosn0  +  ^^yf(^r)sinn0].  (37b) 


The  problem  is  to  minimize  the  functional: 


m  _  m  _ 


(38) 


7=1 


7=1 


The  first  terms  in  each  bracket  are  the  actual  velocities  from  the  elastic  solutions  and  the 
second  terms  are  the  specified  velocities  from  initial  conditions. 


We  have: 


dF 

dAn 


=°'Jr=°- 

dAi2 


dF 

dBn 


=  0,^  =  0, 


dB, 


n 


dF 
dB, , 


=  0 


iZ. 


(39) 


where  i  =  0,L2,...,N. 

The  final  matrix  form  of  Eqn.  (38)  is  obtained  as  follows  by  Least-Squares  procedures. 
Matrix  formulations  are  shown  in  Appendix  C: 


PrP  +  QrQ 

PrC  +  QrD 

"A" 

Pr  QT~ 

'S' 

CTP  +  DTQ 

CtC  +  DtD 

B 

CT  Dr 

T 

where 


A)i 

B0l 

A)2 

Bq2 

A >/ 

=  A 

Bq  1 

4i 

and 

Bn , 

Ai2 

Bnl 

_Anl  _ 

A. 

(41) 
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X' 

'V 

=  S  and 

A. 

_v 

A, oi 

Pi, 01 

Pun 

•••  Put  • 

”  A,nl  ”• 

Pui 

A, oi 

Pl.Ol 

Pl.U 

P2.11 

Pl,n\ 

Pl,nl 

p  = 

JPm,  01 

Pm,0l 

Pm, 11 

-  Pm,U  • 

Pm,n  1 

P  m,nl  _ 

A,oi 

#1,0  1 

#1,11 

#1,11 

A,nl 

Q\,nl 

A, 01 

#2,0/ 

#2,11 

#2,11 

Anl 

$2,nl 

Q  = 

• 

Jim,  01 

#m,01 

#m,ll 

#m,U 

Qm,n\ 

An,n/_ 

Matrix  C  and  D  have  the  same  pattern  as  P  and  Q. 
The  components  of  the  matrix  P,  Q,  C  and  D  are: 


(42) 


(43a) 


(43b) 


Pjjk  =yO)ikUl3\aikrj)smWj 

(44a) 

4jjk  =  y6)lkV\i)(aikrj)cosiQ] 

(44b) 

(44c) 

(44d) 

1.4  Fourier  Transform  Method 

Another  available  method  to  solve  the  problems  other  than  using  the  Least  square 
procedure  is  the  Fourier  transform  technique.  Since  Hankel  functions  are  complex  valued 
there  is  no  point  in  insisting  upon  the  real  form  of  Fourier  series.  Thus  we  write  for  the 
Fourier  transform  of  the  displacement  and  stress  fields: 
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(45) 


u,  =  -XKW(3l(a'-)  +  Bn(CJ)U?(pr)]eme 

V 

/  n=—°° 

ue  =  -^[/A  (6))y®(«r)  -  («)yf(/3r)]em0 

r 

'  n=— oo 


r*  n--oo 

a,  =  2f  £[. \(a)K(ar)  +  Bn(ro)r’0r)V’ 

T  n--oo 

(46) 

are  =  ^  X[iA„(ffl)7lf  (ar) -  iBn(co)T% (pr)Vie 

r  “ 

where  Uf\ar),  ....  ,  are  the  same  functions  defined  by  Eqns.  (19)  and  (23)  with  a 

positive  sign  only,  in  which  =  7T'1  is  used  for  C'0.  The  conditions  An ( co )  and  Bn{co) 
are  to  be  determined  from  the  boundary  conditions. 

1.4.1.  Boundary  Conditions 

Surface  velocities  are  assumed  to  be  given  on  the  surface  of  the  cavity  r  =  a;  i.e.. 


ur(a,d,t)=Vr(6,t)  (47a) 

ue(a,d,t)  =  V9(d,t).  (47b) 

Expanding  the  Fourier  transforms  of  these  into  Fourier  series  we  have: 

Vr(6,  co)  =  Tn(o))ein6  (48a) 

Vg(e,0))=^Sn(0))eine  (48b) 

where  Tn(co)  and  Sn(co)  are  determined  through: 

(49) 

U«»)J  2x{\v,(9,a»y 


where  n  =  0,±1,±2,... 
Since  velocity  becomes: 
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(50) 


Vr(d,a)  =  - Xh'®A  {co)U?\ota) - icoBn (co)Ui3) (pa)]ein6 


a~ 

1 


Ve(6,(0)  =  -  ^coAn(CD)V\\aa)  -  coBn(co)V?\Pa)]eM , 
a„TL 


the  boundary  conditions  from  Equation  (47)  now  give: 

~U?«xa)  U?\Pa)l(An(CQ)]  =  a_ 
yf\aa)  VfXpa)\\Bn((D)]  V 


iTn(co) 

.•W. 


The  solution  of  this  set  is 


An(co)  =  (-)[-iT.(fi>)Vf(/k)  -  Sn((D)U?\pa)]/ det 

CO 

Bn  (©)  =  {-)[Sn{co)Uf\aa)  -  iTn(co)Vf(aa))/da 
CO 

where 


det  =  -Ul3)(aa)V?\pa)  -  Vf\aa)U(3\pa) 


(51) 


(52a) 

(52b) 


(53) 


1.4.2.  Displacements  and  Stresses 

Substituting  the  above  An(co )  and  Bn(co)  into  Eqns.  (45)  and  (46)  we  obtain  the 
solutions  as  follows.  Displacement  fields  are: 


where 


ur(r,  e,co)  =  -  JyTn(coyycor)  +  S„(co)uf\cor)]eM  (54a) 

Y  n~- «» 

ue(r,  6,(0)  =  -%[Tn(co)vycor)  +  iSn{co)v^\cor)]eM  (54b) 

F  n=-oo 


u^(cor)  =  (^)[-Vi3)Q Qo)U(3\ar)  -  V?  (aa)U{3)  (pr)]  /det 


uf  (cor)  =  (^-)[-U?\Pa)U?\ar)  +  U(3\aa)Uf  (pr)]  /det 


v"\cor)  =  (^)[Vf  (pa)Vl3\ar)~  Vi3\ad)V? (pr)] /det 
v(f>(cor)  =  (^)[-U,3)(pa)Vl3)(ar)  -  Uf\aa)V(3)  (Pr)]  jc let. 


(55) 
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Stress  fields  are: 


or(r,  6,(0)  =  ^-  X[iT„  (6))C(cor)  +  Sn((0)t«\(Or)]eM 

F  n=—°° 

CT»(r,  6,(0)  =  ^  2[<T,(<»W>«ar)  +  S,(ffl)j:”(fflr)k“ 

T  n=-« 

^(r,  0,©)  =  ^7-  +  £S.(©)tf  (fflr)]eM 

^  n— — °® 

where 

C«or)  =  (-)[-VfO 3a)Tlf(ar)  -  V?(aa)1™(Pr)]  /det 

f®(<yr)  =  (-)[-U?(Pa)Tl?(ar)  +  Uf\aa)T\^0r)]  /det 

affcyr)  =  (-)[-Vf03a)r‘f  (or)  -  V?(aa)T!»(Pr)]  /det 

(a*r)  =  (-)[-C/f03a)r®  (ar)  +  ^(aa)7^Cj8r)]  /det 

<‘>(6)/-)  =  (-)[V2(3,03a)r<f  (or)  -  Vf(aa)7^03r)]  /det 

<2,(y)r)  =  (-)[-17®03fl)7?(ar)  -  U?\ad)T%  (fir)]  /det . 
ty  / 

Velocity  Boundary  Conditions 

When  we  have  the  boundary  conditions  as 

(i)  Vr(a,0,O  =  K(0,O  =  Vo  for  -0O  <  0  <  0O 
=  0  for  elsewhere 

and 

(ii)  Vs(a,  0,0  =  0, 

then 

ZK  -ea 


(56a) 

(56b) 

(56c) 


(57) 


(58a) 


(58b) 
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(59a) 


=  —sin  (nd0) 

Mt 

£„(©)  =  0.  (59b) 


We  obtain  the  displacement  and  stress  fields  by  substituting  T  (co)  and  S  ( co )  into  Eqns. 
(54)  and  (56)  as  follows: 


and 


u  (r  Q  co)  =  * f  i V" sinwg"  t a 1  t 

r  ’  ’  r„±L  mt  co  U\i\aa)VfHfia)+Vf){oui)Uf\pa) 

1  ^  V^innd,  a .  -Ff(/fa)F,(3)(ar)  +  V™  (cca)V?  (pr)  , 


ug(r,d,co)  =  -Y 

r  z~L  nit 


( ) 1 — 1 — : - 1 — : = — -£ 

co  Uf\aa)Vf\pa)  +  Vf]  (aa)Uf  (pa) 


(60) 


a r(r,0,G))  =  ^y- jN 


V^nnd,  a  Vf (ftQTff (ar)  +  V^aa)Tl(M 
ffl  Uf\aa)Vf  (fia)  +  Vf\cm)Uf\pa) 


Mt 


ae(r,d,co)  =  ^^i 

•  n=~ 


,  Vn  sinn0„  a  Vf> (gg (or)  +  Vf\cm)T^(Pr)  e 
co  )  C/,(3)  (aa)V?  (pa)  +  V?\aa)U™  (Pa) 


UK 


a  (r  e  co)  2/1  f  V"  sin‘ ”flt  ( a  1  zSWZ£W  +  V?{aa)T»{fir) 
rS'  ’  ’  '  r2  -1  ™  TJ^CmnWA^CR^  J.  \/0‘)frf^\TT0)(R„\ 


Jnd 


mt  co  Uf)(aa)Vf(po)  +  V™ (aa)U? (pa) 


(61) 


where  Uf)(ar),....,T<^(Pr)  are  the  same  functions  defined  by  Eqns.  (19)  and  (23)  with 
positive  signs  only,  in  which  C®  =  //(l)  is  used  for  Cjf.  Now  the  inverse  of  Fourier 
transforms  of  these  equations  gives  the  complete  solutions. 

2.0  Elastic-Plastic  Boundary 

The  main  problem  is  to  determine  the  contour  C  that  separates  the  plastic  from  the 
elastic  region  such  that  the  displacements,  as  well  as  stresses,  are  continuous  throughout 
the  exterior  of  the  cavity  surface.  If  this  elastic-plastic  boundary  C  has  been  found,  the 
problem  divides  into  a  pure  plastic  problem  and  a  pure  elastic  one  so  that  the  elastic  and 
plastic  solutions  can  be  obtained  for  each  time  step. 

We  have  three  cases  as  we  evaluate  the  Von-Mises  yield  condition  along  the  each  sector: 


(Case  1) 

|(crr-tJ0)2+6c7re2<2y2 

:  elastic 

(Case  2) 

|(<Tr-CTe)2+6<7re2=2y2 

:  elastic-plastic  boundary 

(Case  3) 

^(<jr-<je)2  +6crre2>2Y2 

:  plastic. 
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As  we  keep  checking  the  yield  condition  exterior  of  the  cavity  surface,  the  elastic-plastic 
boundary  has  been  found  if  the  solutions  at  the  node  point  become  elastic.  Then  the 
stresses  inside  the  boundary  C  are  plastic  and  can  be  obtained  by  using  the  equilibrium 
equations  and  yield  condition. 


elastic 

region 


Figure  3.  Elastic  and  plastic  regions. 


3.0  Plastic  Region 


3.1  Stresses 

We  have  the  equilibrium  equations,  in  which  the  body  forces  are  neglected,  for  the 
plane  strain  case  in  the  cylindrical  coordinates  as  follows: 


i£fi.  =  0 

dr  r  dd  r 


dOre  1  do, 

dr  r  d6 


+  — <7 


rB 


=  0. 


(62a) 

(62b) 


At  the  elastic-plastic  boundary,  plastic  solutions  are  the  same  as  elastic  solutions  because  of 
the  continuity  condition.  Plastic  solutions  at  the  j  =  s  node  on  contour  Ts  now  can  be 
found  with  the  stresses  at  the  previous  contour  FJ+i  from  the  elastic-plastic  boundary  up  to 
the  cavity  wall. 

( Gr  )/•>■  ~  (  Gr  \  [  1  (Gre)  ,.+■+!  ~  ( ^re  )  j'*l  |  ( ' °r  \ /»■  ~  j  Z?.  Vj  -  Q  (63a) 

Ar  r.„,  AO  r;S+1 


(o re).^-{ore).s  1  Wrl+l-((y9)r,  ,  2  ,  _A 

- 1 - 1  AO  +  1  U 

Ar  r.„,  AO  r  ,+1  J 


(63b) 
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Then  we  have  three  components  of  stresses  as  follows: 


Wf  =  Wf*  +~Je^a^r  4i  -(Gre)r\+j-[(<yr)ri  -K)y,J 


,  .  ,  ,  Ar  r .  .  .  .  .  2 Ar , 

(or9)  =  (°J  >  + - —  C(cra)  ]  + - (c) 

7  7  r.„,A0  7  +1  7  r  7 


(64a) 

(64b) 


(<*•),  =  (<T),,  +A/^y2  ~'4(%)/  • 


(64c) 


3.2  Displacements 


In  the  plastic  region,  the  displacements  can  be  determined  from  the  fact  that  er  +  e9 
is  still  given  by  Hooke’s  law,  since  the  sum  of  the  plastic  part  of  er  and  e0  vanishes. 
Assuming  small  strains,  the  displacements  equations  become: 


and 


du  u  1  dv 

—  H - 1 — 

dr  r  r  dd 


l-2v 

2  G 


(Or  +  °e) 


dv  V  |  1  du  _  Gre 
dr  r  r  dd  G 


(65a) 

(65b) 


Then  we  have: 


V1  ur  ,  Vi  V1  ,  V-  _  i  -  2  v 

I  - - 1 - - 


Ar 


rrAQ 


2  G 


{(crr)  -(cra)  _+1} 


(66a) 


V  .  —  V  U  ,  —  U  , 

_V  _V4i  y 


y.  (y),< 


_  _ 


rs+1A0  r.„, 


(66b) 


Now  we  have  two  components  of  displacements  as  follows: 


1  1 
- u  = - XI  -f- 

Ar  Jl  Ar  7'1+1 


V+, “V  .  y.  l-2v 


2G 


{K),+1-K),J  (67a) 


=J_V  <V+.-V'  y 

Ar  J!  Ar  7"‘  r,+1A0  r„,  /x 


(67b) 
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Appendix  A 


Displacement  in  terms  of  Ank  and  Bnlc 
Aside  from  Qxp(-icot),  since 


®(r,d)  =  '£AnHinl\arynS 
«= 0 


'¥(r,0)  =  '£BnH«)(Pr)eine, 
«= o 


we  have: 


=  t,  k-nHlu«xr )  +  or )]«'"“ 


rt=0 


*=0 


fF=i®-4?w>v 


n=0 


dr 


«=0 


V; 


<90 


=  ^iBnknH^(Pr)eine. 


«= 0 

Then  we  have: 


„r=^+i^ 

<9r  r  00 


=  -IXH^'W)  +  +  /5,„n//r)(i3r)}e'"e 


rc=0 


=  -I[^^.(3)(«0  +  /5„tC/f)()3r)]e 

r  n= 0 

1  <90  dy/ 


in  9 


9  r  dd  dr 


=  -  £  RrfW)  +  (/3r)]}e"ifl 


r  “ 


n=0 


=  -  E  [M*  V<3)(ar)  +  Bnk  V?\Pr)]eine . 


r~ 


n=0 
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Appendix  B 


Stresses  in  terms  of  Ank  and  Bnk 

Plane  strain  stress-strain  relationship  is: 
c 7X  =  2Gex  +  X(£x  +  £y) 

(Jy  —  2  G£y  +  X(£x  +  £y) 

oz=X(£x+£y)=v(<Jx  +  oy) 

Vxy  =  Gr*y 


where  shear  modulus  G  and  the  quantity  X  are  the  Lame  constants. 
In  cylindrical  coordinates. 


_  1  dug 

~^~de 


£z=0 

_  1  .1  dur  |  due  ug 

r6  2rd6  dr  r 

£rz  =0- 

The  stresses  in  terms  of  displacements  ur  and  ue  are: 


<J,  =  (2  G  +  X)^  +  M^+^) 

dr  r  dd  r 

ffe  =  (2  G  +  *)(I^  +  i)+^ 

r  dO  r  dr 

dr  r  dd  r 
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„,l  dur  diia 

Trg=G( — L  +  — s- 

r  dO  dr 


(Tn=0 

a0z=O. 

Using  the  above  equations,  we  can  obtain  the  stresses  in  the  cylindrical  coordinates  in 
terms  of  the  unknown  parameters  \k  and  Bnk.  Aside  from  Qxp(-io)t),  since: 


=  ^I(4[-«rfl!')i(°!r)  +  (n2  +«  -  a2r2)H^\ar)] 
+iBnk[nprH^  (fir)  -  ( n 2  +  n)H^firWn6 


^7  =  ^{i\k[nurH"\(ar)  -  (n2  +  n)H^(ar)} 

+BJfirH%lx(fir)  -  (n2  +  n  -  p2r2)H^fir)]}ein6 

r)u  1  °° 

=  ■ -XRJnarC(ar). -  n2H"\ar)} - Bnkn2  H?  firWne 
vu  r  n=0 

r)u  ^  °° 

=  p-X(-A,,»J«r(ar)  +  iS„,  [-,« (0r)  + 

_  ? 

the  final  stresses  become: 

<?r  =  ^tXX^IX"2  +  n  ~ \p2r2)Hn)(ar)  ~  ctrH„\(ar)\ 

r  n=0  *= 1  ^ 

+iB„[-(n2  +  n.  Hi"  (fin  +  nprH"  (/)r)])e" 

-(n2  +  « -  — /32r2  -  a2r2)H(nl)(ar)  +  arH^(ar)] 

r  n~0  *=1  ^ 

+iB,i  t(n2  +  n)H"'([ir)  -  rc/J/tf  “(WDe" 

[-(n2  +  n)H(nl\ar )  +  narH^\  (ar)] 

r  n=0*=l 

+B„,[-(n2  +  n  -  ^2r2)H^fir)  +  (y3r)]}e"*0 

^  =^X2A,(aV-I/JV)//<l)(a^ 

r  n=0  *=1  ^ 

(Tfe  =  0. 
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Appendix  C 


Matrix  formulation  for  Least-Squares  procedures 

Define: 


(I) 


dF 


=  0  case: 


(II) 


=  0  case: 

dB, 


v^djAM-v^hj 

V9j(rrepA,,Bik)-Ve=kr 


f(Ml+M2l=o 

ft  3A, 

A  dh.  dk 
ft  3A„  9A„ 


ft  3B,  3B„ 

■£0l.p-+kp-)=o 

ft  ’  9B,  '  SB, 


(1)  For  A.k 
Define: 

dhj 

dK 

dkj 

Then, 

\\  K  K  ...  UP  +  h  *2  *3  ...  fcm]Q  =  [°  0  0  -  °] 


=  ^(r;’6>>) 


where, 
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» 


A, 01  "■  Pi, 01  Pi,  11  P\,\l  A.nl  Pl.nl 

A, oi  *'■  A,o;  A.ii  P2.11  P2.ni  P2.ni 

P  = 

_Pm,  01  Pm.Ol  An, 11  Pm.ll  Pm.nl  Pm.nl 

#1,01  #1,0/  #1,11  ■"  #1,1/  #l,nl  #l,n/ 

#2,01  $2,01  Q2.ll  ■"  #2,1/  ■"  #2,nl  ”*  #2,n/ 

Q  = 

_#m,01  Qm.Ol  Qm.W  Qm.ll  Qm.nl  #m,n/ 

Also  define: 

^  *3  -  =  H 

[*,  fc,  ...  t.]  =  K. 

Then  we  have: 

HP  +  KQ  =  [0  0  0  ...  0]. 

Transposing  gives: 

(HP  +  KQ)r  =  [0  0  0  ...  Of 
(HP)r  +  (KQ)r=[0  0  0  ...  Of. 

We  now  have: 
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Since  we  have: 


A' 

A 

"V 

A‘ 

y* 

v* 

K 

y* 

k3 

— 

Vri 

— 

K 

and 

— 

y* 

= 

A. 

A. 

A. 

A. 

A. 

we  finally  have: 


A' 

i 

A 

i _ 

y* 

vr2 

y* 

A 

r 

+  Qr 

=pr 

+  Qr 

A 

A_ 

i 

•  ^ 

*  J 
_ 1 

(2)  For  Bik 
Define: 


dhj 

3kj 


=  cjArrdP 

=  djJk{rrdj). 


Then, 

A  K  K  -  KjC  +  lk,  k2  k3  ...  *m]D  =  [0  0 


Matrix  C  and  D  are  same  patem  as  P  and  Q .  Then, 


HC  +  KD  =  [0  0  0  ...  0]. 

Transposing  gives: 

(HC  +  KD)r  =  [0  0  0  ...  Of 
(HC)r  +  (KD)r  =  [0  0  0  ...  Of. 

We  now  have: 
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We  finally  have: 


CrHr  +  DrKr 


0 

0 

0 


"V 

~K 

"o' 

A 

A 

0 

A 

+  Dr 

A 

— 

0 

K. 

_K_ 

0 

1 

_ 1 

X' 

v*. 

v? 

+  Dr 

=  cr 

+  Dr 

V H 

A. 

f5  • 

1 _ 

_vrm_ 

^1  • 

a 

1 - 

If  we  define  the  column  vector  of  unknown  constants  Aik  and  Bik  as 
respectively,  i.e.. 


An 

A>. 

Aa 

Ai 

Ai 

=  A 

A  i 

4,. 

and 

Bnl 

4,2 

Bn  2 

.At. 

Bni. 

we  have  the  velocity  matrices  expressed  in  terms  of  A  and  B  as  shown  below: 


and  B, 
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Now  let's  also  define  the  column  vector  of  the  specified  velocities  at  each  node  point  as 


Then  two  final  matrix  forms  are: 

(1)  For  Art  : 

Pr(PA  +  CB)  +  Qr(QA  +  DB)  =  PTS  +  QrT 
(PrP  +  QrQ)A  +  (PrC  +  QrD)B  =  PrS  +  QrT 

(2)  For  4  : 

Cr  (PA  +  CB)  +  Dr(QA  +  DB)  =  CrS  +  DrT 
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(CrP  +  DrQ)A  +  (CrC  +  DrD)B  =  CrS  +  DrT 


For  both  equations,  the  matrix  form  becomes: 


PrP  +  QrQ 

PrC  +  QrD 

A' 

Pr  Qr 

‘S' 

CtP  +  DtQ 

CTC  +  DTD 

B 

CT  DT 

T 
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